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PREFACE  i 

"Underway  on  nuclear  power."  With  this  simple  but  historically 
significant  message,  the  U.S.S.  Nautilus  heralded  a  new  age  in  power 
generation.  Powered  by  the  first  mobile  nuclear  reactor,  she  set  the 
stage  for  dramatic  developments  in  military  and  industrial  fields.  A 
powerful  vessel,  but  unwieldy  and  undesirably  large  for  a  submarine, 
her  size  was  largely  determined  by  the  size  of  the  power  plant  com- 
ponents. One  of  the  more  significant  of  these  components  is  the  pres- 
surizer  for  the  reactor  primary  coolant  system.  This  tank  is  extremely 
large  due  to  the  conservative  design  necessitated  oy  ignorance  of  the 
thermodynamic  transient  behavior  in  its  pressure  range. 

The  objective  of  this  thesis  is  to  produce  design  data  for  steam 
pressurizing  systems,  Dy  electronic  analog  simulation  of  the  thermo- 
dynamic transients  which  occur  in  the  pressurizer  vessel. 

The  writers  wish  to  thank  Professor  Eugene  E.  Drucker,  and  Profes- 
sor Hugo  U.   Martinez  of  the  U.S.  Naval  Postgraduate  School  for  their 
assistance  and  enthusiastic  cooperation  in  the  prosecution  of  this 
project.         • 
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CHAPTER  1 
INTRODUCTION 

The  prototype  (Lark  1)  of  the  Submarine  Thermal  Reactor  (STR)  was 
the  first  full  scale  power  reactor  in  the  world  to  be  completed  and 
successl'ully  operated.  The  first  central  station  nuclear  power  plant 
in  the  United  States  will  be  the  Pressurized  Water  Reactor  (P'.VR), 
nearing  completion  at  Shippingport,  Pennsylvania.  The  many  highly 
desirable  aspects  of  the  pressurized  light  water  reactor  cause  it  to  be 
one  of  the  most  promising  types,  to  date. 

A  basic  requirement  for  these  light  water  reactors  is  that  a  very 
high  pressure  is  maintained  on  the  primary  coolant.  Pressures  in  the 
neighborhood  of  2000  psia  and  higher  are  presently  in  use.  High  pres- 
sures permit  the  use  of  sub-cooled  water  at  high  temperatures  in  the 
reactor  without  danger  of  boiling. 

When  any  control  program  other  than  constant  average  primary  loop 
temperature  is  used,  a  change  in  the  volume  of  the  primary  coolant  is 
to  be  expected  for  a  corresponding  change  in  power.  A  volumetric  tran- 
sient is  consequently  induced  and  a  surge  tank  is  required  in  the 
system.  Since  the  primary  coolant  must  be  maintained  at  a  high  pressure 
the  surge  vessel  is  also  used  to  perform  this  function. 

The  pressurizer  must  be  strong  enough  to  maintain  steam  and  water 
in  equilibrium  at  high  pressures.  It  must  be  large  enough  to  absorb 
normal  and  accidental  surges  without  permitting  excessive  pressures 
in  the  primary  loop.  Although  thermally  insulated,  it  must  have  internal 
heaters,  capable  of  maintaining  the  water  and  steam  at  the  saturation 


temperature.  These  heaters  must  have  additional  heating  capacity  cap- 
able of  generating  steam  at  sufficient  rate  so  as  to  prevent  excessive 
pressure  drop  in  the  primary  coolant  during  negative  (out)  surges. 

Specifically,  in  this  thesis  we  will  assune  a  primary  loop  working 
pressure  of  2000  psia,  normal  volume  surges  of  4  cubic  feet  and  accidental 
surges  of  7  cubic  feet  as  representative  values.  At  steady  state  water 
and  steam  are  maintained  at  saturation  temperature  (636  °F)  in  the  tank, 
but  only  the  pressure  (2000  psia)  is  transmitted  to  the  cooler  primary 
loop  via  a  relatively  long  stand-pipe.  Assuming  that  the  installed 
heaters  adequately  compensate  for  negative  surges,  this  type  surge  is 
not  considered  a  design  limitation.  Therefore  only  those  surges  caused 
by  an  increase  in  primary  loop  volume,  or  positive  surges,  will  be 
considered  here. 

When  designing  or  sizing  a  pressurizer,  limits  must  be  imposed  on 
the  various  properties  of  the  fluids  in  the  tank.  The  maximum  allowable 
pressure  change  must  be  determined  from  consideration  of  the  strength 
of  the  members  involved.  The  expected  change  in  average  density  and 
therefore  volume  of  the  primary  coolant  must  be  computed.  The  duration 
of  the  surge  is  determined  from  changes  in  power,  either  accidental  or 
deliberate.  The  one  remaining  variable,  the  size  of  the  tank  must  now 
be  determined  on  the  basis  of  the  aforementioned  parameters. 

The  pressure  change  caused  by  a  positive  surge  may  be  minimized  by 
introducing  a  portion  of  the  surge  water  into  the  top  of  the  tank  as  a 
finely  divided  spray.  The  most  conservative  basis  for  pressurizer 
design  is  to  assume  that  no  fepray  is  employed,  the  steam  is  dry  saturated, 
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and  undergoes  isentropic  compression.  Should  the  steam  not  be  ary,  the 
resulting  final  pressure  would  be  lower  than  in  the  case  of  dry  steam. 
The  optimum  and  therefore  minimum  size  required  can  only  be  positively 
determined  after  a  complete  understanding  of  the  thermodynamic  and  heat 
transfer  mechanisms  involved. 

In  this  project,  it  is  felt  that  valuable  information  is  available, 
both  in  the  fields  of  design  limitation  and  in  determining  the  charact- 
eristics of  high  pressure  saturated  steam.  A  Boeing  Electronic  Analog 
Computer  with  associated  4-channel  Sanborn  Recorder  (see  Fig.  1)  is 
utilized  to  simulate  the  pressurizer,  with  the  above  objectives  in 
view. 

An  important  assumption  which  is  made  at  the  outset  is  that  the 
driving  function,  specific  volume  (which  is  directly  related  to  tank 
level)  follows  a  cosine  curve.  Analysis  of  actual  tank  transients 
indicates  that  this  approximation  is  a  good  one  for  a  large  portion  of 
them.  This  type  of  function  can  be  easily  produced  on  the  analog  com- 
puter, while  production  of  more  exact  functions  is  difficult. 

In  the  analysis  of  results,  every  attempt  is  made  to  utilize  the 
principle  of  geometric  similarity,  to  more  nearly  generalize  these 
results. 


BOEING  ELECTRONIC  ANALOG  COMPUTER  AND  SANBORN  RECORD  a 

Figure  1 
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CHAPTER  11 
DEVELOPMENT  OF  THEORETICAL  EQUATII 

The  initial  step  in  the  analog  simulation  oi'  this  problem  consists 
of  a  description  of  the  system  in  terms  of  theoretical  equations. 
.Vhenever  possible  known  thermodynamic  relationships  are  used,  and  when 
this  is  not  possible  empirical  relationships  are  developed.  No  reference 
can  be  found  in  the  literature  treating  the  following  development  of 
descriptive  eolations ,  and  it  is  believed  that  this  approach  may  be  a 
new  one. 

As  an  initial  simplifying  assumption,  an  equivalent  system  is  devised 
which  lends  itself  more  readily  to  analysis  than  does  the  actual  tank. 


\v 


HEAT   SINK  \V 

HIM  It  I 


PERFECT 
INSULATION 


W*~> J 


Equivalent  System 
Figure  2 


The  actual  tank  is  well  insulated,  therefore  the  concept  of  perfect 
insulation  for  the  equivalent  arrangement  does  not  introduce  significant 
error.  The  thermodynamic  system  is  defined  as  the  mass  of  steam  in  the 
tank.  In  the  actual  tank  tnere  is  installed  internally  a  spray  and 
degasifier  assembly.  In  this  analysis  of  positive  surges,  without  spray, 
its  only  effect  is  to  function  as  a  heat  sink.  The  heat  sink  in  the 
equivalent  arrangement  represents  the  tank  walls,  spray  and  degasifier 
assembly,  the  mass  of  water  in  the  tank,  and  all  other  apparatus  within 
the  tank  insulation.  The  piston  represents  the  water  surface,  which 
in  the  actual  tank  performs' work  upon  the  steam  during  an  in-surge. 
Thermodynamic  equilibrium  is  assumed  at  all  times. 

Equation  #1  is  the  well-known  work  equation: 

w  s  -If  dv 

where     w  -  work  (BTU/lbj,  positive  when  work  is  done  on  the  system 
P  =  system  pressure  (psia) 
v  =  specific  volume  of  steam  (cubic  feet/lb) 

Considering  units, 

w  =  -144 |P  dv  -     -.1853JP  dv 

In  Heaviside  operator  form, 

w  -  -if.  1853  P(dv)dtl-  -.1853 JP  (dv)  dt 
pL      dt   J       J   dt 

Equation  #2  is  the  First  Law  of  Thermodynamics: 

Ae  =  q  +  w 
or        e  -     q  ■+-  w  +  e 
where     e  -  specific  internal  energy  (BTU/lb) 

q  =  heat  flow  (BTU/lb),  positive  when  into  the  system 


Equation  #3  is  the  equation  of  state  for  a  non-ideal  gas: 

Pv  =  c1  Z  T 

where  Z  is  defined  as  Pv     (dimensionless  compressibility  factor,) 

C]T 
T   -     absolute  temperature   (°A) 

3 
cl  z  1545  s     .5957  lbf  ft  (gaa  constant) 

13.01x144  iri     lb.°E 

ill 

Equation  #4  relates  P  to  e  and  v  empirically  from  steam  table 
values  (see  Appendix.  1): 

P  =  8.51  e  +  52,590  v2  -  30,700  v  -  3161 

Equation  #5  relates  Z  to  P  and  v  empirically,  in  the  saturated 
and  superheated  region,  1'rom  steam  table  values(see  Appendix  1): 
Z  a  1.715  x  10~4  P  -  9.47  v2  +  6. 043  v  -  .5693 
Note:  Equation  #5  is  not  required  nat hematic ally,  but  is  desirable 
for  computer  circuit  simplicity. 

Equation  #6  relates  heat  flow  to  temperature  differences  and  is 
developed  as  follows: 

First,  consider  the  sink  as  a  system, 

d^sink  =  msink  Csink  ^sink 
where     '^sink  =  Quantity  of  heat  energy  (BTU) 

msirk  =  mass  °f  sifik  (lb) 

Tsink  s  temperature  of  sink  (°R) 

csink  =  sink  heat  caPacity  (BTU/lb°a) 
Rearranging , 

dTSink  -  =Arj    dJ3ink 
"HainkCsink 


Integrating, 

Tsink  =    . — I -isink      +     (T0^si. 

msinkcsink 

but  Qstea7i  *     "  '^ink  =  Q 

and      ^steam  ■  (To)sink  "     To 
Then, 

Tsink  =  _^ Q  +  T0 

msink  csink 

I'Jow,  considering  the  original  system  (the  steam): 

d^  .  -K  (T  -  Tsink) 
dt 
K  represents  the  effective  thermal  conductance  of  the  system  boundary, 

and  will  be  considered  a  constant  in  the  equivalent  system  for  any  given 

surge.  The  units  of  K  are  (BTU/sec°Ii)  •  K  is  ordinarily  the  product 

hA,  when  each  is  determinable. 

Substituting  for  T  .  .  , 

d£  =  -  K  (T  +     ^         -  T0) 
dt  asinkcsink 

Since  Q  =  m  q  ,  where  ms  =  mass  of  steam  (lb) 

da  s   -  K   (T  -  Tc)  -      K       q 
dt      ms  msir.kcsink 

The  problem  which  is  now  encountered  is  that  of  evaluating  the 
equation  constants.  Since  accurate  design  data  are  not  available,  these 
constants  must  necessarily  be  evaluated  in  an  arbitrary  and  approximate 
manner. 

For  the  subject  tank,  m  .  ,  can  be  taken  as  3000  lb.  This  is 
approximated  as  follows: 


mass  of  metal  ~   2000  lb. 

mass  of  water  «   1000  lb. 

Total  3000  lb. 

The  heat  capacity  of  water  at  the  pressure  and  temperature  range 

involved  is  approximately  2.0  BTU/lb°R  .  The  heat  capacity  of  stainless 

steel  (tank  wall  material)  is  approximately  .13  BTU/lb  R  .  A  weighted 

average  heat  capacity  of  the  sink  is  computed  as  follows, 

c  .  .  =  2  x  1000   +   .13  x  2000   .  .75  3TU/lb°R 
slnK  3000 

The  mass  of  steam  in  the  system  will  depend  on  the  initial  tank 
level,  and  saturation  conditions,  but  rarely  varies  greatly  from  an 
average  value  of  175  lb.  Properly,  the  mass  of  steam  should  be  computed 
and  Equation  #6  modified  for  every  run  simulated,  but  for  simplicity 
this  average  value  is  used  for  all  surges.  In  all  simulated  surges, 
rasink»  csink»  anc*  ms  are  taken  as  constants,  while  the  quantity  K  is 
adjusted  to  cause  computed  transient  properties  to  coincide  with  experi- 
mental data. 

Equation  #6  with  constants  evaluated  becomes, 

da,  -  -  5.72  x  10"3  K  (T  -  T  )  -  .445  x  10~3  K  q 
dt 

In  Heaviside  operator  notation,  the  equation  takes  the  form, 

q   =  -  I  C 5-72  x  10"3k  (t  -  V  +  -^  *  10"3k  0 

? 
During  all  runs  the  assumption  is  made  that  the  steam  mass  remains 

constant  throughout  the  surge.  Since  the  energy  introduced  as  work 

during  a  surge  is  on  the  order  of  3000  3TU,  the  maximum  mass  change 

would  be  in  the  neighborhood  of  6  lb.,  if  all  of  the  energy  were  involved 
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in  a  phase  change.  Therefore,  the  percentage  mass  change  can  never  be 

a  significant  value,  especially  since  heat  flows  abundantly  to  the  sink. 

Any  phase  change  will  significantly  affect  K  however,  because  of  the 

heat  transfer  mechanism  involved  (heat  of  vaporization). 

The  six  theoretical  equations  as  developed  and  in  their  most  useful 

forms  appear  as  follows: 

#1       w  =  -  1  j  c12  P  (dv)  dt  j 
p  L.     at    J 

,/2       e  =  q  ■*  w  +  e0 
#3       T  t  P  v 

Ct,  Z 

2 
#4      P  =  c2  e  -  Co  v  +  c^  v   -  c5 

#5       Z  .  c^  P  t  c„  v  -  cg  v   "  c9 


#6       q  =  -  1  Pc   K  (T  -  TQ)  +  c±1   K  q  "j 
where, 


P  L 


c1  z   .5957  c?  =  6.048 

c2  ±  8.51  cg  -  9.47 

c0  r  3.07  x  104  cn  =  .5693 


3  "  "  '  9 

3  -  5.72  x  10 
-3 


c  =  5.259  x  104  c1Q  =  5.72  x  10"- 


cr  =  3161  cxl  =  .445  x  10 

c6  =  1.715  x  10"4  c   =  .1853 
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CHAPTi^t  III 
APPLICATION  OF  EQUATIONS  TO  THE  ANALOG  COMPUTER 

After  the  derivation  of  the  six  theoretical  equations,  a  computer 
circuit  must  be  designed  to  accomplish  their  simultaneous  solution.  In 
the  circuit  design  many  things  must  be  considered;  such  as,  scaling 
factors,  follow-up  time  of  multiplier  components,  amplifier  drift,  volt- 
age magnitudes,  etc.  A  certain  amount  of  latitude  is  allowed  in  circuit 
component  selection,  but  definite  limits  are  present  due  to  stability 
considerations. 

The  six  basic  equations  in  computer  notation  appear  as  below: 


1  [l.lll  (.IF  dv  dt  )  ! 
p  1-         dt     J 


ffl  w  c 

!f2  e*  -  .0333  q   +   -0333  w   +   e 

#3  T  «  3.3574(~j 


ifk  P  -     4.255  e  -  5.1167  v  +■   4.3825  (.02  v^j  -  52.6833 


lb  Z  =  .5145  e  +  3.024v   -  2.3675  (.02  v2)  -  2S.465 

#6       q  s  -  lf.1716  KT  -  .1716  K  TQ  +   .000445  Kql 
P 

The  bar  over  a  quantity  represents  the  quantity  as  a  voltage.  The 
detailed  conversion  of  the  theoretical  equations  to  the  above  "machine" 
equations  comprises  Appendix  II. 

Fig.  3  shows  the  resulting  machine  circuit  in  block  diagram  form. 
For  the  detailed  circuit,  see  Appendix  II.  Briefly,  the  circuit  is 
designed  to  function  as  follows: 
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In  the  driving  function  circuit  a  cosine  function  is  added  elect- 
rically to  a  constant  voltage,  the  resultant  voltage  transient  very 
closely  approximates  the  specific  volume  throughout  a  surge.  The 
amplitude  of  the  cosine  function  represents  one-half  the  change  in  v  , 
while  the  period  of  the  function  is  twice  the  time  of  the  surge  (a  t) . 

Within  this  same  circuit  the  time  derivative  of  specific  volume,  dv 

dt 
is  also  generated. 

A  function  multiplier  is  used  to  obtain  the  quantity  lr t   while  a 


servo-multiplier  produces  the  products  P  dv  and  P  v.  An  integrating 

dt 

circuit  integrates  P  dv  as  a  function  of  time  to  produce  specific  work, 

_  dt 

w  .  Another  integrating  circuit  produces  the  time  integrated  specific 

heat  flow,  q  ,  using  temperature  inputs. 

A  summing  amplifier  sums  e0,  w  ,  and  q  to  produce  e.  Another  sums 

e,  v,  and  v^  to  produce  P.  A  third  sums  P,  v,  and  v  to  produce  Z. 

A  dividing  circuit  performs  the  division  P  y  to  produce  T~. 

Z 
The  circuit  required  to  solve  the  set  of  equations  consists  of 

eleven  amplifiers,  four  sign-changing  amplifiers,  two  function  multipliers, 
one  duo-channel  servo-multiplier,  and  twenty-one  50,000  ohm  potentio- 
meters. The  net  result  is  a  rather  complex  circuit,  but  one  which  does 
not  seem  to  be  subject  to  further  simplification  without  sacrificing 
necessary  accuracy.  Since  instability  in  this  type  q£   computer  is  a  direct 
function  of  the  number  of  components  used,  circuit  simplicity  has  been 
a  constant  objective. 

This  phase,  the  designing  and  setting  up  of  the  circuit,  requires 
a  great  deal  of  computer  experience  and  proficiency.  For  the  uninitiated, 
many  trial  and  error  situations  are  encountered,  and  with  a  circuit  of 
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this  complexity  the  process  can  be  very  time  consuming. 

The  following  chapter  will  deal  with  the  problem  of  producing 
quantitive  results,  assuming  the  circuit  has  been  designed,  balanced 
and  tested. 
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CHAPTLR  IV 
ANALOG  COIJPUTEil  SOLUTION  OF  THEOKhriCAL  EQUATIONS 

After  the  analog  computer  circuit  is  designed,  assembled,  and  rough 
qualitative  results  obtained,  the  next  phase  consists  of  producing 
accurate  quantitative  solutions.  The  first  step  in  accomplishing  this, 
consists  of  static  adjustment.  With  various  initial  values  of  e  ,  and 
v  set,  the  e,  P,  Z  and  T  computing  circuits  are  adjusted  to  give  steam 
table  values  for  these  quantities  over  the  entire  anticipated  range. 
This  step  of the  co  puter  set-up  is  very  exacting  and  time  consuming 
but  is  a  necessary  step  prior  to  adjusting  for  dynamic  accuracy. 

Prior  to  the  dynamic  adjustment,  a  hand  solution  of  the  system  of 
theoretical  equations  is  obtained  for  use  as  an  adjustment  reference. 
A  set  of  initial  conditions,  and  a  driving  function  are  arbitrarily 
selected,  and  using  the  Heun  Methodof  numerical  integration,  a  point- 
by-point  solution  of  the  equations  is  calculated  (see  Appendix  III). 
The  various  thermodynamic  quantities  as  calculated  are  plotted,  the 
resulting  curves  serving  as  the  final  reference  for  dynamic  adjustment 
(see  Fig.  4) • 

Next,  the  computer  is  made  to  duplicate  the  hand  solution.  This 
is  done  by  adjusting  the  w  and  q  integrating  circuits  until,  the  w  and  q 
curves,  as  viewed  on  a  properly  calibrated  Sanborn  Recorder,  are  of  the 
proper  shape  and  magnitude  (see  Fig.  5).  These  adjustments  in  general 
are  minor. 

The  value  of  careful  static  adjustment  is  now  apparent  since  no 
further  adjustment  of  the  e,  P,  Z  and  T  circuits  is  found  necessary. 

1   Special  case  of  Runge-Kutta  Method;  see  Kopal,  Numerical  Analysis, 
pp.  202-205,  Wiley  (1955). 
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The  curves  of  these  quantities  are  1,01V  found  to  match  the  correspono 

hand  solution  curves  very  well  (see  Figs.  4  and  5/. 

Once  the  hand  solution  is  duplicated,  a  reference  for  measuring  K 

is  now  available.   In  the  hand  solution,  the  arbitrarily  selected  value 

of  K  is  .9254.   The  corresponding  setting  of  potentiometer  a,,  is  now 

J-4 

readily  observed  (usually  in  the  neighborhood  of  .100,  varying  slightly 
from  day  to  dayj .  iiince  K  and  the  a.,  setting  are  proportional,  K  can 
be  determined  at  any  time  by  the  simple  relation, 

*=  KoxCajJ 

It  is  now  possible  to  vary  K  over  a  large  range  of  values  by  merely 
adjusting  a  single  potentiometer,  thus  varying  the  integrated  q  curve, 
and  all  other  curves  of  thermodynamic  properties.  By  simulating  the 
driving  function  for  a  given  surge,  K  can  be  varied  until  the  proper 
pressure  curve  is  obtained  (i.e.  a  pressure  curve  which  matches  the 
pressure  curve  resulting  from  an  actual  tank  transient;.  The  ne^t 
chapter  discusses  the  investigation  of  actual  surges  by  simulator,  to 
determine  values  of  K  for  these  surges.  Correlation  of  K  with  surge 
types  is  also  discussed. 
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CHAPTER  V 
EXPERIMENTAL  RESULTS  AND  CONCLUSIONS 

Appendix  IV  is  a  compilation  of  data  obtained  from  the  actual  surge 
bank.  Eighteen  positive  surges  are  represented,  spray  being  utilized 
in  only  one.   It  is  evident  that  possession  of  many  more  runs  would  be 
desirable,  but  the  present  number  must  sui'fice  siace  ...ore  could  not  be 
obtained.   Fortunately  these  eighteen  surges  represent  a  1'air  cross- 
section  of  positive  surges  encountered. 

The  procedure  for  simulating  a  given  surge  is  outlined  in  Chapter 
IV  ana  briefly  is  as  i'ollowo: 

A  proper  driving  function  (specific  volume  J  is  set  up  on  the  comput- 
er, providing  for  the  generation  of  a  cosine  curve  of  proper  amplitude 
and  period.  Initial  values  of  v,  e,  T,  and  a  trial  value  of  K  are  set, 
and  a  surge  is  generated.  K  is  adjusted  for  several  trials  until  the 
generated  property  curves  match  the  corresponding  actual  tank  curves 
(see  Fig.  6  for  Sanborn  recording  of  this  procedure). 

Each  of  the  eighteen  surges  is  simulated  in  this  manner,  and  the 
proper  value  of  Redetermined  (see  Table  1).  kany  attempts  at  correlat- 
ing K  v/ith  various  quantities  (a  v,  av  ,  av_  ,  etc.)  were  fruitless,  but 

vo  AXj 
a  fair  correlation  with  surge  duration  {/ntj   is  possible  (see  Fig.  7). 

It  is  likely,  in  fact,  that  K  is  a  function  of  several  or  all  of  the 

variables  considerea,  but  the  effect  of  each  seems  small  in  comparison 

to  the  effect  of  At.   3y  simple  curve  fitting,  a  correlating  equation 

is  derived, 

K  =  3-75  *  25 

t  -  17 
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Nothing  i3  claimed  for  this  equation  except  that  it  describes  a  K  which 
aids  in  predicting  pressure  rises  on  the  safe  side  (i.e.  in  most  casus 
the  predicted  pressure  rise  will  be  somewhat  larger  than  the  actual 
pressure  rise).  Column  4  of  Table  1  lists  the  K  for  each  surge  as 
predicted  by  formula. 

In  general,  it  can  be  observed  that  K  does  not  vary  greatly  for 
any  surge  of  over  30  seconds  duration,  and  that  it  becomes  nearly  cons- 
tant for  surges  longer  than  90  seconds,  it  is  also  true  that  pressure 
rise  is  not  sensitive  to  moderate  changes  in  K.  The  last  column  of 
Table  1  lists  errors  encountered  when  a  constant  value  of  3-75,  the 
limiting  value  of  K,  is  used  in  the  computer. 

The  mechanism  of  heat  transfer  for  the  actual  tank  can  only  be 
surmised.   However,  it  is  logical  to  expect  greater  turbulence  for  short- 
er, faster  surges,  resulting  in  somewhat  larger  values  for  K.  This 
could  be  due  to  a  larger  heat  transfer  surface  between  the  steam  and 
water  in  the  tank.  The  effectiveness  of  spray  in  increasing  K  is  readily 
observed  in  run  1506,  Table  1.  Whereas  a  K  of  approximately  5-0  could 
have  been  expected  without  spray,  the  actual  K  observed  is  13.8  . 

From  the  general  characteristics  of  the  K  vs.  Surge  Duration  curve 
(Fig.  7;,  it  appears  that  K  approaches  a  limiting  value  of  about  3.75 
for  surges  of  long  duration.  For  geometrically  similar  tanks  utilizing 
the  sane  pressure  range,  the  same  limiting  value  of  K  could  reasonably 
be  expected. 

Realizing  that  a  correlation  exists  between  K  and  the  driving 
function,  it  is  now  possible  to  predict  pressure  rises  for  all  conceiv- 
able driving  functions.  This  is  done  in  the  following  manner: 
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as  m  : 1  i  ^ 


SUM  NO. 

3U      b 

(sec; 

K  (exp.; 

K  (for  . , 

2045 

75 

4.31 

4.17 

7 

2302 

120 

4.12 

4.00 

-  4 

1630 

120 

3.96 

4.00 

-12 

2123 

100 

4.26 

4.05 

-  5 

1653 

95 

4.26 

4.06 

18 

2229 

90 

^•07 

4.07 

-22 

2023 

65 

2.73 

4.27 

-34 

1556 

65 

4.54 

4.27 

12 

2010 

60 

3.15 

4.32 

-21 

1101 

27.5 

6.30 

6.15 

30 

2240 

45 

6.25 

4.65 

12 

2255 

32.5 

5.55 

5.35 

18 

2248 

35 

3.90 

5.15 

7 

1911 

37.5 

4.36 

5.00 

66 

2113 

37.5 

6.50 

5.00 

60 

2039 

45 

5.3^ 

4.65 

66 

2247 

45 

6.02 

4.65 

90 

1506 
(spray) 

45 

13.3 

4.65 

200 

■   i_,rror  in  predicting  oressure  rise  (in  psia)  when  using  a  constant 
value  of  K,(  the  limiting  value  of  3.75;,  rather  than  the  exper- 
imental or  formula  values. 
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A  driving  function  is  e:  tablished,  the  proper  value  of  K  is  set, 
and  the  resulting  pressure  rise  recorded.  In  an  attempt  to  use  the 

'inciple  of  geometric  similarity,  diraensionless  parameters  are  used 
when  possible.  Accordingly,  a  family  of  cui'ves  is  obtained  relal 

a  v  ,  aj_  f  and  at  (see  Fig.  8 ) . 

v     ?n 
o        .    o 

since  the  actual  tan1;  surges  vary  fro:j  27.5  secor.ds  to  120  seconds, 
the  predictions  ax-e  only  valid  within  this  range.  However  it  appears 
that  safe  extrapolations  may  be  made.  An  attempt  to  justify  liraitin 
values  of  the  prediction  curves  is  based  on  the  following  reasoning: 

Initially  the  system  consists  of  dry  saturated  steam.  As  a  surge 
progresses  work  is  done  on  the  system,  and  heat  will  flow  out  provided 
a  temperature  difference  exists.  Examination  of  Fig.  9  shows  that  the 
isentropic  compression  process  prescribes  the  practical  upper  limit  of 
over-pressure.  It  is  believed  that  a  process  following  the  saturation 
line  prescribes  the  lower  limit  of  over-pressure  for  this  tank,  oince 
the  surge  introduces  colder  water  into  the  tank,  a  potential  sink  exists 
which  could  conceivably  cause  the  process  to  enter  the  wet  steam  region. 
however,  a  careful  study  of  actual  surges  indicates  that  the  stea.i 
remains  superheated  throughout  an  in-surge.  Very  long  surges  appear  to 
rallel  closely  the  saturated  line  with  slight  superheat.  In  addition, 
all  surges  appear  to  be  initially  isentropic  in  nature . 

Based  on  the  aoove  reasoning  ana  actual  tesc  results,  it  is  believed 
t  the  prediction  curves  approach  finite  values  as  1 1  approaches  in- 

_  dty.   rhe  limiting  value  of  a_t   (for  any  given  A_v  )  is  that  which 

P?  vo 

would  be  obtained  in  a  process  following  the  saturation  line. 
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TEMPERATURE  -ENTROPY     CHART  (APPROXIMATE) 

FIGURE   9 
26 


For  very  short  surges,  the  Limiting  condition  appears  to  be  i 

of  isentropic  compression,  thus  establishin  a 

for  any  given  Av  .  On  the  actual  installation,  surges  of  less  than  10 

v 
seconds  duration  <.re  unlikely  uuc  to  primary  loo.)  circulation  ti.ie 

considerations.  Thus  the  extrapolated  curves  for  values  of  a t  from  10 

seconds  to  30  seconds  would  seem  reasonable.  Using  the  prediction 

curves  (Fig.  S) ,  Table  2  shows  a  comparison  of  predicted  and  actual 

pressure  rises. 

It  is  believed  that  basic  objectives  of  this  project  have  been 
realized;  however  many  avenues  for  further  study  appear  to  be  open. 
Although  a  procedure  for  rational  analysis  has  been  indicated,  it  appears 
that  a  study  of  the  effect  of  spray,  and  of  negative  surges, by  these 
methods  would  be  of  value. 

In  the  field  of  heat  transfer,  a  great  deal  remains  to  be  done. 
The  combining  of  all  heat  transfer  mechanisms  into  a  single  coefficient 
is  perhaps  an  over-simplification,  although  effective  in  this  analysis. 
It  is  known  that  several  thermal  paths  exist,  and  a  more  complete 
analysis  could  perhaps  isolate  the  effects  of  each.  Heat  can  flow  from 
steam  to  metal,  from  steam  to  \vater  (directly  or  by  condensation;,  and 
from  metal  to  water.  A  comprehensive  simulator  study  coula  provide 
knowledge  in  this  relatively  unexplored  field. 
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TA3LJ  2 


RUN  NO. 

2302 

2045 
2123 
1556 
2229 
1506 

2023 
2010 
2243 
2240 
2247 
2113 
1101 
2255 
1911 
1658 
1630 
2039 


(;>sia) 


179 
317 
310 

236 

298 
302 
214 
212 
257 
130 
380 
150 
139 
179 
420 
335 
329 
362 


ACTUAL  I 

(psia; 

174 
313 
300 
229 
303 
239 
242 
232 
261 
110 
314 
187 
112 
152 
417 
310 
324 
341 


PREDICTION 

__(j2sia_i 


5 

4 

10 

7 

-  5 
13 

-28 
-20 

-  4 
20 
66 

-37 

27 

27 

3 

25 

5 

21 
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APPENDIX  1 

DERIVATION  cF  EMPIRICAL  RELATIONS  FOR  PRESSURE 

AND  COMPRESSIBILITY  FACTOR 

Since  e  can  be  computed  from  the  First  Law,   e  ■  q  +  w  +  e 
and  v  is  available  as  a  driving  function,  it  is  convenient  to  express 
P  as  a  function  of  e  and  v  .  This  can  be  done  by  constructing  a 
empirical  relationship  from  steam  table  values,  which  will  oe  accurate 
within  certain  limits.  The  pressure  range  which  is  considered  signifi- 
cant for  this  problem  is  1700  -  2700  psia,  and  no  accuracy  is  claimed 
beyond  these  limits. 

Since  e  cannot  be  taken  directly  from  the  st-eam  tables  in  the 
superheat  region,  a  plot  of  enthalpy  vs.  pressure  is  made  initially, 
showing  lines  of  constant  temperature  and  constant  superheat  (see  Fig. 
10). 

Next,  at  various  temperatures,  various  values  of  h  ,  P  ,  and  v 

are  listed.  For  each  of  these  values  e  is  computed  from  the  relation 

e  =  h  -  144  P  v  (see  Table  3)«  From  this  taole,  a  new  chart  is 

J 
constructed  relating  e  ,  P  ,  v  and  T  (see  Fig.  11).  Using  this  chart 

as  a  basis,  initially  an  expression  for  e  =  f  (P,v)  is  obtained  in 

the  following  manner: 

The  average  slope  of  constant  v  lines  is  .1176.  For  any  one  line, 

e  =  .1176  P  +  c-i  1    .   The  constant  c-y  is  determined  as  a  function  of 

v  (see  Table  4  and  Fig.  12;,  c-^  -  1562  v  +  533.  The  resulting 

expression  for  internal  energy,  e  »  .1176  P  +■  1562  v  f  538  ,  when 

tested  at  the  extreme  limits  is  found  to  jive  considerable  error. 

Although  undesirable  from  the  standpoint  of  computer  circuitry,  a.. 
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TABLii  3 


r  (°f) 

630 

(sat. ) 
640 

(sat.) 
650 

(sol. ) 
660 

(sat. ) 
670 

(sat.) 

650 

(sat. ) 
690 

700 


h   (3TU/lb) 


1156.1 
1167.7 
1145.8 
1141.1 

1185.7 
1167.0 
1145.6 

1130.6 

1167.0 

1146.3 

1121. U 
1118 . 5 

1167.7 
1147.3 
1123.8 
1104.4 

1150.0 
1123.2 
1099 . 5 
1087.7 

1132.3 
1107.1 

1072.5 

1067.2 

1156.6 
1137.3 

1114.3 

1160.6 
1142.5 


P  (psia) 


1700 

1900 

191 V. 3 

18U0 
1900 
2000 

2059.7 

2000 
2100 
2200 
2203.2 

'2100 
2200 
2300 
2365.4 

2300 
2400 
2500 
2531.8 

2500 
26U0 
2750 
2708.1 

2500 
2600 
2700 

2600 
2700 


v  (ft3/lb) 


.2553 
.2296 

.5540 
.1992 

.24U7 
.2168 
.1936 

.1798 

.2053 

.1546 

.1633 

.1616 

.1962 
.1763 
.1575 
.1442 

.1702 
.1526 

.1342 
.1277 

.1484 
.1319 
.1137 
.1115 

.1594 
.1447 
.1299 

.1549 
.1415 


P  v 


8u.5 
76.6 
71.8 
70.3 

30., j 
76.4 
71.7 
68.7 

76.2 
71.8 

06.4 
66.2 

76.3 
72.2 

67.3 

63. 0 

72.5 

67.5 
62.2 
59.7 

68.7 
63.5 

36.5 
56.0 

73.3 
69.3 
64.8 

74.6 
70.7 


e  (aru/ib) 


1105.6 
1091.1 

iu74-o 
1070.3 

1105.4 
1090.6 
1073.9 

1061.3 

1090.8 
167^.5 
1054.6 
1052.3 

1091.4 
1075.6 

101; 
1041.4 

1077.5 
1060.4 

1037 . 6 
102  . 

1063 . 6 
1043.6 
1016.6 
1011.2 

1082.8 
1067.5 
1049.5 

1036. u 
1071.8 
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ui 
(C 

v3 


i    u. 


TaJlK  k 


V 

e 

P 

.1176  P 

C14 

.12 

1021.0 

2620 

308.5 

712.5 

.13 

1030.4 

2510 

295.5 

735.0 

.14 

1038.2 

2409 

284.0 

754.2 

.15 

1045.4 

2312 

272.  G 

773.4 

.16 

1051.5 

2217 

261.0 

790.5 

.17 

1056.8 

2140 

252.0 

804.3 

.18 

1062 .0 

2057 

242.0 

820.0 

.19 

1066.5 

1982 

233.5 

333.0 

.20 

1070.5 

1915 

225.5 

845.0 
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additional  modification  is  necessary,  making  e  =  g(  P,  v  ,  v  ) . 

o 

To  accomolish  this,  the  expression  P  -  S.51  e   ■  c.  ,.  v  ■*■  d  /  v 

15       io 

+  c-^y  is  set  up.  Substituting  actual  steam  table  values  at  three 
points  (  P  &  1900,  2100,  2300)  three  simultaneous  equations  in  cjc  , 
C-j/  ,  and  c,r,  result.  Solving  for  these  constants  and  substituting, 
the  final  equation  for  P  is  : 

P  -  8.51  e  +  52,590  v2  -  30,7^0  v  -  3161 


The  calculation  of  an  empirical  relationship,  l.   =.  f  (p,  v)  follows 
similar  lines.  First  Z  vs.  P  is  plotted  from  the  steam  tables, 

showing  lines  of  constant  v  (see  Fig.  13) .  The  average  slope  of  the 

4 
constant  v  lines  is  1.715  x  10""  ,  the  equation  for  any  one  line  being 

Z  -  1.715  x  10~4  P  +   clc,  .    c-ja  is  determined  as  a  function  of 

v  ,  c. d  =  2.943  v  -  .318  (see  Table  5  and  Fig.  14) .  The  expression 
18 

for  Z  as  a  function  of  P  and  v  is,   Z  =  1.715  x  10   P  +  2.943  v 
-   .318 

Once  again  a  v^  term  is  found  necessary"  for  desired  accuracy,  and 
using  the  same  procedure  as  before  the  resulting  expression  for  Z  is: 
Z  =  1.715  x  10"4  P  -  9.47  v2  +      6.048  v  -  .5693 
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TA3LE  5 


-4 
1.715  x  10   P 


'18 


.12 

.13 
.14 
.15 

.16 

.17 
.18 
.19 
.20 
.21 


.4640 
.4845 
.5030 
.5225 
.5380 
.5510 
.5640 
.5770 
.5890 
.6005 


2613 
2509 
2410 
2312 
2211 
2140 
2060 
1985 
1915 
1850 


.4475 
.430u 
.4130 
.3960 
.3790 
.3670 
.3530 
.3405 
.3280 
.3170 


.0165 
.0545 
.0900 
.1265 
.1590 
.1840 
.2110 
.2365 
.2610 
.2835 
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APPENDIX  11 
DERIVATION  OF  MACHINE  EQUATIONS 


Scaling  factors  are  assigned  as  follows: 


max 


3000 

50 


60 


cK 


_       V 


:nax 


raax 


50 


.01 


*  p    =    T« 


aax       r     1500 


o< 


e     - 


CX. 


Z      - 


o< 


W 


c*\ 


oT 


T 
max 

50 

e:oax 

e 

max 

— 

1^00 
50 

^raax 

^uax 

= 

1 

50 

vv 
'  'max 

s 

%L 

w 
max 

50 

;jmax 
%iax 

= 

£0 

50 

T       = 
t 

1 
1 

= 

30 


3o 


,02 


=        1 


Conversion  of  Theoretical  Equations  to  Machine  Equations 


#1 


w    =     -  1 

P 


c12  f 


'(f) 


dt 


W       x 


-  1    1  .1351  (.1  P  dv)  dt]   10  c*pc*v c*t 
p     L  dt  '      J      o/  ^ 


°C.<*t 


W       s       - 


1.111     (.1  F 
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dv  \  dtl 
dt  /      J 


#2 


e     -     q    ¥    w    +    e( 


e     = 


2a  q 


+     ^ 
c*0 


w       +  cx-e 


cxTe      e0 


=     .0333     q      +     -0333    w 


#3 


P  v 


.01)T 
.02/ 


o^n    o^         .01   P    V 

ci^z^r  1.02  z 


.01  F   v  1 
.02  L       J 


#4 


P     -s       c2  e     -     c^  v    +  c^  v 


-     c, 


P     =  o<_  c0     e 


<*, 


-<*v  c< 


C*r 


+    50  c^y^v  ( .02  v^)      -  < 


o<; 


3 


P     =  4.255     e 


-  5.1167  v     +     4.3825  (.02  v2)-  52.6833 


#5 


Z     =     c6  P    +-  c?  v    -     c8 


c^c^p     P  -»-  c7cXv 


V2 


-  50  cfio^v  c*v  (.02  v  J    -  cs. 
— rr o<7. 


o< 


Z      = 


.5145  P    +    3.024    v 


xz 


-     2.3675   (.02  v   )     -  28.465 


#6 


-  lfci0K   (T  -  T0)+    cu  K     qj 


q    =  - 


lphoKoCp   (rng  _+.  c,lk 


*q°<t 


<*t 
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-1 


q  =  -  1  T.1716  K  T   -  .1716  K  Tn  +   .0)0445  K  q 
p  L 


Calculated  potentiometer  settings,  and  components  seiecied: 


al 

=  D.  F. 

depei 

ident 

a8  = 

.4255 

a2 

11 

n 

aQ  = 

.6864 

a3 

11 

■ 

n 

aiu  = 

.5117 

a4 

=  .1111 

aii  ■ 

.6333 

a5 

r  .3333 

a12  * 

.3U24 

a6 

=  .3357 

ai3  = 

.5145 

■7 

=  .3333 

a14  = 

K 

R   =  1  M 

R,  r     m     10  Li 

45 

He  -  i0  M 
af   =   1  m 

Roo  z  .1  M 
Rf30  s  lu  M 
R6o  ..ill 

R5   a  .1  M 

R^j_  =  10  M 
R68  -  .1  M 


cfl  = 

,!/£ 

cf43- 

2/£ 

Cf8  = 

,1/Af 

CflO« 

,l/*f 

Cf23- 

.ou5/f 

a15  = 

.4383 

ai6  = 

.6864 

a!7  = 

.5000 

ai8  = 

.0089 

a19  = 
a20  = 

L 

45 
.2368 

a21  = 

.ii82 
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R-IOM 
-AVvWv        *■ 


R?IM 


■  io^=  ioU.sin  m^ 


R  =  /M 

-www- 


Q:  -.ll 


■lou.u)  sin  cut: 


-10  u.  cos  wt 


R<MM 
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Rf-IM 


-10 


R„-2M 
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-.ipjr 


at 
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R-IM 

-vwim - 
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APPEI'JDU  111 
HAND  SOLUTIOM  OF  THEORETICAL  EQUATIONS 

The  first  step  in  obtaining  a  hand  solution  consists  of  combining 
the  six  theoretical  equations  to  form  two  dependent  differential 
equations : 

Substituting  from  equation  #2  into  equation  #4,  and  then  into 
equation  ,;/'l  one  obtains, 

$  =  -clA[v^  +  e„-  |v  +|  v*.  fVjdy 

Since  \j    =       U  +      u,     Co$    cot 

then  -^    -      -U0Co   S/ilwt 

gl«    ulU)S(;iatWtt^w.^  +  &  ./  +  e„  -  &  "I 

Substituting  from  equations  ,74  and  //5  into  equation  ff3 ,  the 
resulting  expression  for  T  becomes, 

it  is  now  convenient  to  assi-r:  values  to  the  driving  function, 
initial  conditions,   and  to  evaluate  the  constants: 
uQ     =     .019  ft3/lb  c2     =     8.51 

OJ      =     .0418  rad/sec  C3     -     3.U7  x  llA 


vmean  =     'l72  rt'ib  c<+     =     5.259  *  L°U 


eo 


1066  3TU/lb  cr  =  3161 

cx  =  .5957  c6  =  1.714  x  10~4 


CY 

as 

6.U48 

c8 

s 

9.47 

cy 

= 

.5b93 

c12 

^ 

.1853 
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Evaluating, 

i* 


=    /.253  *l0*%iy\  .o4/8 "t  [j2  +"U7  -  36«0  V    +  £>i6G  v*  -fr  6  9*] 


Equation  #6  is, 

^    r    -CWKT    +    C,0KT0  -  C„K^ 

Arbitrarily  assigning  a  value  to  K  ,  evaluating  constants,  and 
substituting  the  already  derived  expression  for  T  into  equation  //6, 
K  =  .9254 

clu  =  5'72  X  i0~3 
c±i  =  4.4i>  x  10~4 

T0  =  10940  °a 

^  -  -5/.8v  fa  -fxa  -36jo  v/  fb!8ovz+$>94]  _<?- 

a<  k_$__  d-  4-/2mo  q  tS.770 


\A    r  W   +  540  V  -  3/0  i/2-  t  30+  j 


The  two  resulting  dependent  first  order  differential  equations 
bee one , 

where  /C,  ^  /C2  ,  /f3  /^  are  tabulated  values   of      V       ,   ^he  driving  function 
(see  Table  C) . 
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TABLr,  6 


Time 

UQ  COS  cut 

V 

*L 

K2 

K3 

H 

0 

.0190 

.1910 

-9.9o 

230.5 

395.7 

1    0 

5 

.U186 

,ly06 

-9.88 

231.  U 

395.5 

2.6l  x 

10"^ 

10 

.0174 

.1893 

-V.S2 

232.0 

395.1 

5.10  ' 

15 

.ol>4 

.1873 

-9.71 

235.  o 

394.3 

7.33  ' 

20 

.0127 

.1847 

-9. 58 

239.0 

393.1 

9.33 

25 

.0095 

.1815 

-9.41 

242.5 

391.8 

10.87  ' 

30 

.U059 

.1779 

-9.24 

247.5 

3V0.3 

11.90  ' 

35 

.0020 

.1740 

-V.02 

253.  o 

388.6 

12.47  ' 

•  40 

-.0020 

.1700 

2^v. 5 

330.6 

12.47  ' 

45 

-.OU59 

.1561 

-8.62 

264.5 

385.2 

ll.9o  ! 

50 

-.u095 

.1615 

-8.37 

272.3 

383.1 

1U.B7   ' 

55 

-.0127 

.1593 

-8.26 

276.0 

382.1 

9.33  ' 

60 

-.0154 

.1567 

-8.13 

2al.O 

3&1.0 

7.38  • 

65 

-.0174 

.1547 

-o.02 

283. 7 

380.1 

5.10  ■ 

70 

-.0186 

.1534 

-7.V5 

285.0 

379.6 

2.61  ■ 

75 

-.olyO 

.1530 

-7 .  94 

236.5 

379.5 

0 

80 

-.ui86' 

.1534 

-7.95 

28^. Q 

379.6 

-2.61  ' 

e5 

-.0174 

.1547 

-8.02 

233.7 

380.1 

-5.10  ' 

90 

-.0154 

.1567 

-8.13 

28i.O 

381.0 

-7.33 

95 

-.0127 

.1593 

-8.26 

276.0 

382.1 

-9.33  ' 

100 

-.0095 

.1615 

-8.37 

272.3 

333.1 

-10.37  ' 
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The  next  step  consists  of  numerically  integrating  the  above 
equations  to  obtain  point -by-,  joint  values  of  a  ,  4SL  ,    vJ    ,  and  d°? 

o    '  US  '  d-t 

The  Heun  L'ethod  of  numerical  integration  is  used,  the  general  form 


bein^: 


where 


VJ     -    a   (  OJ  ,  a  ,  V ) 

%■  --  fa  -f  h  I  fa) 

tOc     --   &>..,  ■+  £  (  ca,„,  +  wc       ) 

where  r  rM  _  ,      ,    ,       . 

h  is  the  interval  between  points  and  is  selected  as  5  seconds. 
Table  7  lists  the  resulting  values  for  Q.  ,  C?x    ,  UJ  ,    and  VJ 
obtained  over  a  time  of  1U0  seconds  (21  points). 

Next,  the  resulting  values  of  Q     and  U  are  substituted  back  into 
the  original  theoretical  equations,  resulting  in  point-by-point  values 
for  the  various  properties  throughout  the  run.   These  values  are  tab- 
ulated (see  Table  8),  and  are  plotted  as  a  function  of  tine  in  Fig.  ks 
Chapter  IV. 
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TA3LZ  7 
HAND  SOLUTION   (K  =.9265 J 


w 


w. 


0 

0 

1 

5 

2 

10 

3 

15 

4 

20 

5 

25 

6 

30 

7 

35 

8 

40 

9 

45 

10 

50 

11 

55 

12 

60 

13 

65 

14 

70 

15 

75 

16 

30 

17 

85 

IS 

90 

19 

95 

20 

100 

0 

-.0273 

- . 1083 

-.3410 

-.8431 

-1.5479 

-2.4138 

-3.4896 

-4.7612 

-6.1402 

-7.5359 

-9.1078 

-10.6548 

-12.1476 

-13 . 4633 

-14.6571 

-15.6849 

-16.5174 

-17.1731 

-17.5063 

-17.4614 


0 

•.OilO 
-.0215 
-.0716 
-.1292 
-.1527 
-.1957 
-.2327 
-.2760 
-.2756 
-.3026 
■.3062 
-.3126 

•.2345 
-.2433 
-.2317 
-.1794 
-.1536 
-.1087 
■.0246 
.0426 


0 

.1509 

.6011 

1.3431 

2.3614 

3.6260 

5.0948 

-.7173 

3. 4311 

10.1536 

11.8070 

13.3039 

14.5556 

15.4969 

16.0765 

16. 2714 

16.0843 

15.5423 

14.6944 
13.6080 

12.3475 


0 

.0604 

.1197 

.1772 

.2302 

.2757 

.3119 

.3372 

.3483 

.3407 

.3206 

.2781 

.2226 

.1539 

.0779 

0 

■.0743 

-.1420 

-.1972 

-.2374 

•.2659 
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TABLE  8 
PROPERTY  VALUES  FRuLI  HAND  SOLUTION 

I 
P       i      7. 


0 

5 

10 

15 
20 

25 
30 
35 
40 
45 
50 

55 
60 
65 
70 
75 
SO 

85 

90 

95 
100 


1066.0 

1066.12 

1066.49 

1067.0 

1067.52 

1068.08 

1068.68 

106V.23 

1069.67 

1070.01 

1070.22 

1070.20 

1069.90 

1069.35 

1068.61 

1067.61 

1066.40 

1065.02 

1063.2V 

1062.10 

1060.69 


.1910 
.1906 
.1893 
.1873 
.1847 
.1815 
.1779 
.1740 
.1700 
.1661 
.1615 
.1593 
.1567 
.1547 
.1534 
.1530 
.1534 
.1547 
.1567 
.1593 
.1615 


1966.9 
1969.7 
1986.4 
2015.4 
2047.0 

2086.9 
2134.2 
2190.1 
2243.0 
2297.3 
2361.5 
2391.9 
2427.1 
2447.0 
2459.6 
2458.1 
2440.8 
2390.1 
2370.9 
2323.0 
2282.1 


.5774 
.5774 
.5773 
.5766 
.5760 
.5747 
.5734 
.5718 
.5699 
.5679 
.5653 
.5638 
.5616 
.5597 
.5578 
.5560 
.5546 
.5499 
.5520 
.5520 
.5517 


1092.1 
1091.3 
1093.4 
1098.9 
1101.9 
1106.5 
1111.4 
1118.8 
1123.1 
1127.9 
1132.6 

1134.3 
1137.0 
1135.4 
1135.4 
1135.5 
1133.2 
1128.6 
1129.9 
1125.5 
1121.6 
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APPENDIX  IV 

ACTUAL  TAIJK  TRANSILfoT  DATA 


Run  #2247 

mo 

.  169.06  lb. 

Time 

Pressure 

Level 

Volume 

Spec.  Vol. 

0 

2002 

32.2 

31.7 

.1875 

5 

2012 

32.4 

31.6 

.1867 

10 

2042 

33.2 

31.2 

.1843 

15 

2094 

34.6 

30.4 

.1797 

20 

2168 

36.7 

29.1 

.1720 

25 

2233 

39.7 

27.3 

.1612 

30 

2280 

41.8 

26.1 

.1542 

35 

2305 

43.1 

25.4 

.1500 

40 

2316 

43.7 

25.1 

.1485 

Run  #2240 

0 

5 

10 

15 
20 

25 

30 
35 
40 


m  s  167.46  lb, 


2002 
2007 
2026 
2047 
2071 
2091 
2104 
2112 
2112 


32.8 

33.0 


33.5 


34.1 
34.9 
35.6 

36.3 
36.8 
37.0 


31.4 
31.3 
31.0 
30.6 
30.2 
29.7 
2v.5 
29.0 
28.9 


.1875 
.1870 
.1850 

.1827 
.1803 
.1773 
.1760 
.1731 
.1726 
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Run  #2229 

mo  " 

167.13  lb. 

Time 

Pressure 

Level 

Volume 

Spec. Vol. 

0 

1992 

32.5 

31.5 

.1885 

10 

2007 

33.0 

31.3 

.1873 

20 

2033 

33.9 

30.8 

.1843 

30 

2070 

35.2 

3C 1 

.1800 

40 

2117 

36.3 

29.0 

.1735 

50 

2170 

38.6 

23.0 

.1675 

60 

2220 

40.4 

27.0 

.1615 

70 

2264 

41.8 

26.1 

.1560 

80 

2295 

42.7 

25.7 

.1537 

Run  #2128 

0 

2000 

10 

2012 

20 

2037 

30 

2072 

40 

2114 

50 

2162 

60 

2211 

70 

2244 

80 

2288 

90 

2299 

100 

2299 

m„  =  169.05  lb, 

32.3 
32.7 
33.4 
34.5 
36.0 
37.8 
39.5 
41.0 
42.0 
42.6 
42.7 


31.7 
31.5 
31.0 

3U.4 
29.5 
28.5 
27.4 
26.6 
26.0 
25.7 
25.6 


.1875 
.1861 
.1832 
.1796 
.1743 
.1684 
.1620 
.1572 
.1537 
.1520 

.1514 
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Run  #2113 

mo  = 

172.7  lb. 

Time 

Pressure 

Level 

Volume 

Spec. Vol. 

0 

2000 

31.2 

32.4 

.1878 

2.5 

2002 

31.4 

32.3 

.1870 

5 

2019 

31.8 

32.0 

.1854 

7.5 

205/* 

32.4 

31.6 

.1830 

10 

2090 

33.3 

31.0 

.1795 

12.5 

2124 

34.3 

3^.5 

.1765 

15 

2155 

35.3 

29.9 

.1730 

17.5 

2172 

36.2 

2V.4 

.1703 

20 

2182 

36.9 

29.0 

.1630 

22.5 

2186 

37.3 

28.7 

.1662 

25 

2187 

37.6 

28.5 

.1650 

27.5 

2187 

37.9 

28.4 

.1645 

30 

2187 

38.0 

28.3 

.1640 

Run  #1101 

0 

2.5 
5 

7.5 
10 

12.5 
15 

17.5 
20 


m  =  169.5  lb, 
o 


1990  31.6 

1995  31.8 

2009  32.3 

2039  32.8 

2066  33.5 

2085  3A.4 

2097  34.9 

2102  35.3 

2102  35.5 


32.1 
32.0 
31.7 
31.4 
31.0 

30.4 
30.1 

29.9 
29.8 


.1892 
.1888 
.1870 
.1854 
.1830 
.1793 
.1775 
.1765 
.1760 
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Run  #2255 
Time 

0 

2.5 

5 

7.5 
10 
12.5 
15 

17.5 
20 

22.5 
25 

27.5 
30 

Run  #1506 

0 

5 
10 
15 
20 

25 
30 
35 
40 
45 


m0  = 

:  175.0  lb. 

Pressure 

Level 

Volu  ..ie 

Jpec .Vol . 

1995 

30.0 

33.0 

.1385 

2002 

30.1 

33.0 

.1885 

2011 

30.3 

32.9 

.1880 

2032 

30.6 

32.7 

.1870 

2043 

31.1 

32.4 

.1850 

2055 

31.7 

32.0 

.1828 

2094 

32.5 

^1  A 

31.6 

.1805 

2121 

33.3 

31.1 

.1775 

2134 

34.0 

30.7 

.1755 

2143 

34.6 

30.3 

.1730 

2146 

35.0 

30.0 

.1715 

2147 

35.3 

29.9 

.1710 

2147 

35.5 

29.8 

.1703 

m0  s  167.7  lb. 

1952  30.8 

1958  30.9 

1938  31.7 

2032  33.0 

2082  34.4 

2130  35.7 

2174  37.1 

2217  38.5 

2238  39.7 

2241  40.0 

53 


32.6 
32.5 
32.0 
31.3 
30.4 
29.7 
28.8 
28.0 
27.3 
27.1 


.1944 
.1940 
.1908 
.1867 
.1312 
.1770 
.1713 
.1670 
.1629 
.1616 


Run  #1556         mo  = 

:  159.7  lb. 

Time 

Pressure 

Level 

Volume 

Spec. Vol. 

0 

1393 

31.2 

32.4 

.2030 

5 

1394 

31.3 

32.3 

.2022 

10 

1923 

32.0 

31.9 

.2000 

15 

I960 

33.0 

31.3 

.1962 

20 

1998 

34.2 

3C5 

.1912 

25 

2034 

35.4 

2v.8 

.1868 

30 

2063 

36.5 

29.3 

.1837 

35 

2086 

37.3 

28.7 

.1798 

40 

2101 

38.0 

28.3 

.1774 

45 

2114 

33.7 

28.1 

.1761 

50 

2122 

39.1 

27.5 

.1723 

55 

2122 

39.4 

27.4 

.1718 

Run  #201C 

) 

0 

1930 

5 

1933 

10 

1942 

15 

1962 

20 

1988 

25 

2021 

30 

2057 

35 

2094 

40 

2124 

tf 

2146 

50 

2159 

55 

2162 

m 


o     - 


164.6     lb. 

31.0 

31.0 

31.1 

31.5 

32.3 

33.2 

34.4 

35.4 

36.4 

37.0 

37.6 

37.9 
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32.5 

.1975 

32.5 

.1975 

32.5 

.1975 

32.1 

.1951 

31.7 

.1927 

31.2 

.1398 

30.5 

.1855 

29.8 

.1312 

29.3 

.1782 

28.9 

.1757 

28.6 

.1740 

28.4 

.1727 

Run  #2022 

mo 

-  166.6  lb. 

Time 

Pressure 

Level 

Volume 

opec. Vol. 

0 

1946 

31.0 

32.5 

.1952 

5 

1949 

31.0 

32.5 

.1952 

10 

1962 

31.2 

32.3 

.1940 

15 

1982 

31.7 

32.1 

.1923 

20 

2003 

32.5 

31.6 

.1899 

25 

2040 

33.4 

31.0 

.1862 

30 

2074 

34.5 

30.4 

.1827 

35 

2111 

35.7 

29.7 

.1734 

40 

2144 

36.5 

29.2 

.1753 

45 

2168 

37.2 

28.8 

.1730 

50 

2130 

37.7 

28.5 

.1711 

55 

2186 

33.1 

28.3 

.1700 

60 

2188 

38.3 

28.2 

.1694 

Run  fr'l63C 

)            m 

-  198.3  lb. 

0 

1938 

20.0 

39.0 

.1967 

10 

1952 

20.6 

38.6 

.1948 

20 

1982 

21.5 

33.1 

.1920 

30 

2022 

22.8 

37.3 

.1880 

40 

2074 

24.8 

36.1 

.1830 

50 

2132 

26.7 

35.0 

.1775 

60 

2184 

27.3 

34.4 

.1722 

70 

2224 

29.8 

33.2 

.1670 

BO 

2252 

31.4 

32.2 

.1623 

90 

2262 

32.6 

55 

31.5 

.1589 

rtun  #1911 

;o 

-  134.7   lb. 

Time 

Pressure 

Level 

Volume 

0 

1910 

23.2 

37.1 

2.5 

1920 

23.6 

36.8 

5 

1935 

23.8 

36.7 

7.5 

1958 

24.2 

36.6 

10 

1990 

25.4 

35.8 

12.5 

2027 

26.0 

35.4 

15 

2072 

27.5 

34.5 

17.5 

2117 

28.8 

33.8 

20 

2161 

30.0 

33.1 

22.5 

2207 

31.8 

32.0 

25 

2252 

33.0 

31.3 

27.5 

2287 

34.5 

30.4 

30 

2312 

35.7 

29.7 

32.5 

2323 

36.7 

29.1 

35 

2327 

37.3 

28.7 

37.5 

2327 

37.5 

23.6 

opec .Vol. 

.2007 
.2004 
.1993 
.1974 
.1950 

.1914 
.1872 
.1830 

.1788 
.1737 
.1690 
.1645 
.1607 
.1572 
.1553 
.1548 
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Run  ,t:l65o 
Tine 


ra, 


Pressure 


0 

1972 

10 

1972 

20 

1932 

30 

2002 

40 

2045 

50 

2092 

60 

2147 

70 

2200 

80 

2249 

90 

2273 

95 

2282 

134.1  lb. 
Level 


26.2 
26.4 
26.5 
27.4 
28.9 
30.3 
32.5 
34.6 
36.3 
37.7 


7oln::e 


35.3 
35.2 
35.1 
34.5 
33.7 
32.6 
31.6 
30.4 
29.3 
28.5 
28.2 


8pec  .Vol. 


.1918 

.1916 

.1906 

.1877 

.1830 

.1775 

.1713 

.1651 

.1591 

.1544 

.1532 
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Run  #2248 

mo 

.  158.3  lb. 

Time 

Pressure 

Level 

Voiune 

Spec. Vol. 

0 

1903 

31.8 

31.9 

.2017 

2.5 

1903 

31.8 

31.9 

.2017 

5 

1907 

31.8 

31.9 

.2017 

7.5 

1912 

32.0 

31.8 

.2008 

10 

1923 

32.3 

31.7 

.2000 

12.5 

1948 

32.3 

31.4 

.1982 

15 

1987 

33.6 

31.0 

.1958 

17.5 

2033 

34.6 

30.3 

.1912 

20 

2084 

35.8 

29.6 

.1869 

22.5 

2125 

37.0 

28.9 

.1825 

25 

2150 

38.0 

28.3 

.1787 

27.5 

2162 

38.6 

27.8 

.1755 

30 

2164 

39.2 

27.6 

.1742 

32.5 

2164 

39.4 

27.5 

.1735 

35 

2164 

39.5 

27.4 

.1730 
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Run  #2302 

!            m 

0 

s   162.5  lb 

i 

Time 

Pressure 

Level 

Volume 

Spec. Vol. 

0 

1883 

3-.0 

33.0 

.2032 

10 

1891 

30.2 

32.9 

.2025 

20 

1902 

30.5 

32.8 

.2020 

30 

1918 

31.0 

32.5 

.2000  . 

40 

1935 

31.6 

32.1 

.1975 

50 

1955 

32.3 

31.7 

.1,52 

60 

1975 

33.1 

31.2 

.1921 

70 

1996 

34.0 

30.7 

.1890 

00 

2017 

34.3 

30.2 

.I860 

90 

2038 

35.5 

29.8 

.1835 

100 

2054 

36.1 

29.4 

.1810 

110 

2062 

36.4 

29.3 

.1804 

120 

2062 

36.5 

29.2 

.1798 
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Hun  #20 45 

m0 

-  167.0  lb. 

Time 

■■—  i  ■  ■  i  1 1   i.   .   | 

Pressure 

— — > 

Level 

Volume 

Spec. Vol. 

1 

0 

1932 

32.0 

31.8 

.1904 

5 

1990 

32.5 

31.6 

.1894 

10 

2005 

32.8 

31.4 

.1880 

15 

2025 

33.5 

31.0 

.1856 

20 

2054 

34.5 

30.4 

.1823 

25 

2091 

35.5 

29.8 

.1785 

30 

2130 

36.6 

29.1 

.1746 

35 

2166 

37.5 

23.6 

.1710 

40 

2197 

38.5 

28.0 

.1674 

45 

2227 

39.5 

27.4 

.1640 

50 

2252 

40.5 

26.9 

.1610 

55 

2275 

43-. 3 

26.5 

.1584 

6U 

2290 

41.9 

26.1 

.1562 

op 

2295 

42.1 

26.0 

.1548 

70 

2295 

42.5 

25.7 

'  .1539 
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Run  ;/2039 
Tims 


o 
Pressure 


0 

2.5 
5 
7.5 

10 

12.5 

1.3 

17.5 

20 

22.5 

25 


27.5 


30 
32.5 

35 
37.5 


1981 
1981 
1986 
1995 
2003 
2025 
204^ 
2062 

2137 
2177 
2217 
2252 
2286 
2312 
2322 


165.5  lb. 
Level 


Voluae 


Jpec.Vcl, 


32.7 

3i.5 

.1905 

32." 

31.4 

.1902 

33.0 

31.3 

.1896 

33.2 

31.2 

.1888 

33.6 

JI. 0 

.1874 

34.1 

30.6 

.1854 

3  4.  7 

3o.3 

.1832 

35.4 

29.8 

.x8J2 

36.1 

29.4 

.1770 

37.8 

28.5 

.1734 

33.2 

28.2 

.1698 

3  j.  2 

27.6 

.1662 

40.3 

27.0 

.1626 

41.2 

26.5 

.1590 

42.4 

2>.8 

.1560 

43.1 

25.4 

.1540 
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Simulation  of  steam 
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